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PARTIALLY  CONSTRAINED  IMPINGING  JETS 

ABSTRACT 

Qualitative  explanations  of  the  motion  of  shaped  charge  liners  have 
been  based  on  impact  of  two  plane  Jets  in  which  the  moving  fluid  is 
surrounded  by  four  stagnant  regions,  all  at  the  same  pressure.  Actually, 
the  motion  is  initiated  by  the  difference  between  the  high  pressure  in 
the  detonation  products  on  one  side  of  the  liner  and  atmospheric  pressure 
on  the  other.  This  report  considers  symmetrical  impact  of  two  Jets,  each 
partially  constrained  on  one  or  both  sides,  in  which  the  stagnant  regions 
%re  not  all  at  the  same  pressure. 
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1.  INTRODUCTION 


Treatises  on  hydrodynamics  generally  contain  discussions  of  two  standard 
examples  of  plane  incompressible  jet  flows;  viz,: 

a.  Efflux  of  liquid  with  high  pressure  at  infinity  from  a  reservoir 
with  straight  walls  into  a  jet  surrounded  by  stagnant  fluid  at  lower  pressure; 

b.  Impact  of  two  jets,  in  which  the  moving  fluid  is  surrounded  by  four 
regions  of  stagnant  fluid,  all  at  the  same  pressure. 

Contemplation  of  these  examples  suggests  the  problem,  to  determine  the  fol¬ 
lowing  flow: 

c.  Impact  of  two  jets,  each  partially  constrained  on  one  side  by  straight 

walls,  in  which  the  jets  are  also  partially  bounded  by  four  regions  of  stagnant 

fluid,  not  all  at  the  same  pressure. 

The  symmetrical  form  of  (c),  shown  schematically  in  Fig.  1  (with  an  inflection 
point  A^  on  the  high  pressure  boundary),  will  be  discussed  in  this  note. 

A  perfectly  obvious  generalization  of  (c)  is 

d.  Impact  of  two  jets,  each  partially  constrained  on  both  sides  by 

straight  walls,  In  which  the  jets  are  also  bounded  by  four  regions  of  stagnant 

fluid,  not  all  at  the  same  pressure. 

The  symmetrical  form  of  this  flow  is  discussed  at  the  end  of  this  note  by 
straight  forward  modifications  of  the  mathematical  apparatus  used  to  describe  (c). 

The  flows  to  be  constructed  are  interesting  for  their  own  sakes  as  jets 
that  can  be  described  explicitly  In  relatively  simple  terms.  Additional  in¬ 
terest  might  be  stimulated  by  the  following  considerations.  Flows  of  type  (b) 
have  been  used  to  suggest  a  qualitative  explanation  of  the  motion  of  the  liner 
of  a  shaped  charge.  In  reality,  the  motion  is  produced  by  the  difference  be¬ 
tween  the  high  pressure  in  the  detonation  products  on  one  side  of  the  liner  and 
atmospheric  pressure  on  the  other  side.  Furthermore,  at  an  intermediate  stage 
only  part  of  the  liner  has  collapsed  or  begun  to  collapse,  while  the  rest  is 
still  rigid.  In  an  admittedly  imperfect  way  flows  (c)  and  (d)  more  nearly 
approximate  these  features  than  (b). 
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We  note  also  that  various  pressures,  velocities,  anu  widths  are  referred 
to  in  Fig,  1,  and  four  additional  parameters  will  appear  in  the  discussion 
in  the  following  sections.  Let  us  suppose  we  have  experimental  measurements 
of  the  angles  Qy  0^,  the  three  widths  h^,  hy  and  h^.,  and  assume  we  know  p^ 
(which  might  be  atmospheric  pressure)  and  the  density  p  of  the  jets.  Then 
we  could  apply  the  seven  equations  (2.5),  (3*4)  -  (3*6),  (4.4),  (4.6),  and 
(4.8)  to  determine  the  four  relatively  uninteresting  mathematical  parameters 
and  the  three  important  physical  parameters  p^,  U^,  and  U,_,  (which  will  be 
defined  in  equations  2.4  and  2.5)* 

2.  MATHEMATICAL  FORMULATION  OF  PROBLEM 

Plane  irrotational  incompressible  flow  can  be  characterized  by  a  complex 
potential  function 

(2.1)  $(z)  =  0  +  i*. 

Here  $(z)  is  an  analytic  function  of  the  complex  variable  z  =  x  +  iy,  0  is 
the  velocity  potential  function,  the  stream  function,  and 

(2.2)  w  =  u  -  iv  =  d$/dz 

is  the  complex  velocity.  The  pressure  p  within  the  jets  is  determined  by 
Bernoulli's  equation 

(2.3)  p  +  -  p  |  w  |  *~  =  c  :  stant 
Conditions  on  the  Jet  boundaries  are  characterized  by 

(2.4)  P  =  Pa  ,  a  =  1,  5, 
and  thus  |  w  |  assumes  corre spending  constant  values 

(2.5)  I  V  |=  ua,  a  =  1,  5. 

To  seek  $(z)  or  w(z)  directly  in  the  z -plane  is  hopeless,  since  the 
location  of  the  jet  boundaries  characterized  by  (2.4)  or  (2.5)  is  not  known 
a  priori.  The  classical  artifice  for  overcoming  this  difficulty  is  to  invert 
(2.2)  to  determine 

(2.6)  z  =  f (w ) 

where  f  is  an  analytic  function  of  w.  Then  straight -streamlines  (walls,  or 
the  axis  of  symmetry)  have  as  their  images  in  the  w  -  or  hodograph -plane  seg¬ 
ments  of  lines  through  the  origin,  and  free  jet  boundaries  correspond  to  arcs 
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of  circles  (2.5)  with  centers  at  the  origin.  Thus  the  impinging  Jets  of  Fig.  1 

map  onto  the  interior  of  the  region  shown  in  Fig.  2.  The  circular  cuts  A_A. 

11  ,3  4 

and  A^A^  appear  as  a  matter  of  necessity,  while  the  cut  A^AgA^  has  been  in¬ 
troduced  for  later  convenience. 

Boundary  conditions  and  other  relevant  properties  of  f(w)  can  be  determined 
as  follows.  On  the  straight  streamlines  of  Fig.  1  dz  is  parallel  to  W  and  dw 
is  parallel  to  w.  Thus  on  all  straight  segments  shown  in  Fig.  2 

(2.7)  Im  w2  dz/dw  =  Im  w2f  =  0 

On  the  free  Jet  boundaries,  which  are  also  streamlines,  dz  is  again  parallel 
to  w,  and  or  the  circular  arcs  (2.5)  dw  is  parallel  to  iw.  Thus  on  all  cir¬ 
cular  arcs  shown  in  Fig.  2 

(2.8)  Re  w2dw/dz  =  Re  w2f  (w)  =  0 

To  guarantee  finite  non-zero  jet  widths  at  infinity  in  the  z-plane,  f(v;  should 

1  it 

have  logarithmic  singularities  at  A^,A  ,A^A,-.  At  A^A^Aj^A^  and  at  Ag  the 
function  f(w)  should  be  finite,  and  as  a  matter  of  convenience,  arbitrarily 
choose 

(2.9)  f(o)  =  0. 

3.  CONFORMAL  MAPPING  OF  HOLOGRAPH  IMAGE  ONTO  HALF-PLANE 

As  an  aid  to  determining  the  functional  form  of  w  f  (w)  it  will  be  con¬ 
venient  to  map  the  interior  of  the  curve  shown  in  Fig.  2,  slit  along  A<-Ag,  onto 
a  half  plane.  To  do  this,  first  note  that 

(3*1)  W  =  log  (w/t^) 

maps  the  region  in  question  onto  the  region  with  polygonal  boundary  shown  in 
Fig.  3.  If  we  take  account  of  symmetry,  then  by  the  Schwarz -Christoff  el 
formula, 


(3.2) 


w  =  -  a 


[(l^-l)(|2-a^)(|2-a^)]  0,b 


where  1  <  a^  <  a^<  a<_  will,  for  suitable  choices  of  the  positive  parameters 
a,  & y  a^,  and  a<-  yield  the  desired  mapping  onto  the  upper  half  of  the  £  -plane. 
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In  the  calculation  of  (3*2),  use  that  branch  of  the  integrand  that  is  positive 

for  t  >  ac*  For  later  reference  we  also  note  that 
3 


(3.3) 


1  dv 


/„2  2\ 

a  U  -a^) 


1  *  '  L(?2-l)(C2-3)(?2-4)]  °-5 


(5.4) 


(5.5) 


(5.6) 


(5-7) 


1 

a  I  =  a  J 

a  *2  =  a  I 


(a?  -  I2)  d| 


a  i_  =  a 
3 
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f. 


£ 

a  \  =  a  j" 


[(l-|2)(a2-|5)(a2.|2)]  0.5 

(aj^2)  d4 

[(|2-l)(a2-|2)(a^-|2)] 

0.5 

(52  -aj)  d| 

[(E2-l)(52-a2)(a|-|2)] 

0.5 

i  (S2-aj;)  <34 

[(|2-D(|2-a2)(a2-|2)] 

|  0.5 

=  -  ©, 


=  log  ux/u5 


*  V®4 


=  Jt  +  9 


The  constant  a  can  be  evaluated  as  follows.  Since  the  only  singularities 

of  dW/d  {;  are  at  +  1,  +  a  ,  and  +  a  ,  then 
“  —  3  “  5 

-4ai(i  +  i  +  i. )  -  \  (awM)a*  =  f  (dw/d4)d4 

9  ^Cl+C2+C5  *^C 


where  C.^,  C2,  and  are  paths  shown  schematically  in  Fig. 4,  and  C  is  a  circle 

It!-  const  >  a  .  Clearly 
i  5 
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f  (dW/d4)  <34  =  -  a  I  (| _1  +  ...)d|  =  -  2nia 

J  c  J  o 

Since  by  (3-4)  to  (3-7),  -  4aii  +  I?  +  1^)  =  -  4*1,  this  yields 


(3-8)  a  =  2. 


4.  CONSTRUCTION  OF  f(v) 

2  » /  N 

Recall  that  w  f  (w)  is  alternately  real  or  pure  imaginary  on  the  segments 

i 

of  the  real  axis  of  the  £  plane  with  end  points  and  }  for  y  -  2,  3,  5. 
Thus  it  must  have  branch  points  at  these  places,  and  should  contain  a  factor 


<c2  /2<t2  -  a2)‘/2 

where  r,  s,  and  t  are  odd  integers.  Since  f  should  have  logarithmic  sin¬ 
gularities  at  A^,  Ag  and  A^,  5  =  5,  then  df/d£  should  have  simple  poles  at 

the  corresponding  points.  Since  furthermore  f  should  be  finite  at  k0  and 
A  this  suggests  the  form 


(4.1) 


o  t 

w  f  (w)  = 


ea2  - l)0,5 


tit2  -  a2) 


0.5 


(f 


2w„2  2v  3.5 

-  a,)U  -  a.) 


2  2 

The  factor  £  -  in  the  denominator  will  enable  us  to  include  the  case 

in  the  following  discussion.  Now,  by  (4.1)  and  (5*3) 


(4.2)  df  -2£  _ 1 _ 

w  ttt2  -a2)(s2  'a2) 

where  p  >  0  merely  determines  the  geometrical  scale  in  the  z  =  f  (£)  plane 
Since  w(oo)  =  0,  then  by  (2.9)  f(w(ao))  =  f(o)  =  0,  and 
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(4.5) 


(df/d£)d£. 


The  uniqueness  of  our  choice  for  (4.l)  can  be  shown  as  follows.  Let  us 
multiply  the  left  members  of  (4.l)  and  (4.2)  by  an  analytic  function  H(£). 

2  ’/  v 

To  preserve  the  alternation  of  real  and  imaginary  values  of  w  f  (w)  on  the 
real  axis,  H  must  he  real  there.  To  prevent  the  introduction  of  new  branch 
points  and  singularities,  H  must  have  no  singularities  in  the  closed  upper- 
half  plane,  and  the  analytical  continuation  of  H  into  the  lower  half  plane 
is  also  free  of  singularities.  Hence  H  is  constant. 

It  remains  to  show  that  f(w)  has  the  desired  properties  at  A^,  Ag,  A^, 

and  A  .  First  note  that  t  =  0  corresponds  to  w  =  11 .  Thus  by  (3-l)  and  (3*2) 
y  _L 

£  =  (v  -  ux)  g(w  -  ux) 


where  g  is  an  analytic  function  of  w  -  and  g(o)  =^=  0.  Since  df/d£  has  a 
simple  pole  at  (;  =  0,  then  f(w(£))  has  a  logarithmic  singularity  there,  and 
thus  f(w)  also  has  a  logarithmic  singularity  at  w  =  U^.  A  similar  argument 
determines  the  behavior  of  f  at  £  =  *  a^  or  *  a^,  with  the  unimportant  dif¬ 
ference  that,  for  example,  (3*l)  and  (5*2)  imply 


(£  -  a^)0’5  =  (w  -  U^e  5)h(w  -  U?e  5) 

where  h  is  analytic,  h(0)  =£  0,  etc.  Hence  f(w)  has  the  required  logarithmic 
singularities. 

By  (3.2),  in  the  neighborhood  of  infinity 


where  c^  =  1.  Hence  W  =  -2  log  (;  +  m(l/0  where  m  is  an  analytic  function  of 
l/(;.  Then  by  (3 . l) 

W  »  c  "2  n(l/0 

where  n  is  analytic  and  n(0)  =^=  0.  Then  by  (4.2) 
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dfM  -  i  2  '  5k(l/0 

where  k  is  analytic  and  k(o)  =^=  0.  Thus  f(£)  will  be  regular  at  oo* 

The  width  of  the  jet  at  in  the  z -plane  can  be  determined  by  consid¬ 

ering  the  expansions 


(k.k)  =  2(3  rt/u^  a^ 

and  rate  of  mass  flow 

(4.5)  \  =  2  p  7t  p/a|  a^ 

t 

Similarly,  at  or  A^  we  have  widths 

(4.6)  =  P  n/U5  a^  (a^T  -  a^) 

and  rate  of  mass  flow 

(4.7)  ^  =  p  n  p/a^  (a|  -  a^) 

1 

and  at  A^  and  A^  the  total  width 

(4.8)  y  2p  n/U^  (a|  -  ap 
and 

(4.9)  y  =  2  a  n  p/a^  (a|  -  a|) 

As  we  would  expect  from  the  law  of  conservation  of  mass 

"l*  Hj*  2  “5  • 
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Finally,  the  ratio  of  the  rate  of  mass  flow  at  to  that  at  is 
(4.10)  Mj/H-  =  a^/a*  -  1 

In  Fig.  1  the  straight  walls  were  adjacent  to  the  low  pressure  regions. 
Would  it  be  possible  to  place  them  adjacent  to  the  high  pressure  regions? 

If  we  proceed  purely  formally,  we  merely  have  to  replace  the  simple  poles  of 
df/clt;  at  *  a5  simple  poles  at  *  1.  However,  the  following  intuitive 
considerations  show  that  this  process  leads  at  least  to  an  indeterminacy. 
Suppose  a  flow  of  the  desired  type  exists.  For  our  modification  consider 

d  log  (df/dOM  =  -  v  /w  +  .. 

Near  the  end  of  one  wall  the  behavior  of  this  logarithmic  derivative  is  dom- 

•  , 

inated  by  the  branch  point  at  |  =  a,y  But  for  a^  <  |  <  a^,  Iw  /w  <  0.  Thus, 
as  one  would  expect,  the  streamline  leaving  the  wall  at  A^  bends  initially 
toward  the  low  pressure  region,  as  shown  in  Fig.  5.  Now,  without  changing  the 
flow  field,  we  can  extend  the  walls  into  the  stagnant  high  pressure  regions. 
Since  by  appropriate  changes  of  scale  we  can  always  make  the  gap  between  the 
ends  of  the  extended  walls  be  of  unit  length,  this  means  that  the  location  of 
the  point  of  detachment  is  indeterminate. 

5.  PARTIALLY  CHANNELLED  IMPINGING  JETS 

To  produce  flows  with  impinging  jets  that  are  partially  bounded  by 
straight  walls  on  both  sides,  it  will  suffice  to  replace  the  simple  poles  of 
df/d£  at  +  by  simple  poles  at  *  a^,  where  1  <  <  a^.  Now  (4.2)  becomes 

(5-1)  df  _  2£  _ _1 _ 

«  "  v  '  £  <£2  -  .?H£2  -  a2) 

while  (3.2)  remains  unchanged.  The  presence  of  the  additional  parameter 
will  make  it  possible  to  vary  the  location  of  A^  in  Fig.  6,  for  example,  while 
the  locations  of  A 2  and  A^  and  the  directions  of  the  walls  are  held  constant. 
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The  calculation  of  the  various  jet  widths  and  rates  of  mass  flow  are  perfectly 
straightforward  exercises  which  we  shall  not  repeat.  It  should  be  remarked 
that  if  there  is  an  inflection  point  the  wall  can  be  extended  into  the 
high  pressure  region  again,  just  as  in  the  discussion  of  Fig.  5. 
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FIG.  I  FLOW  IN  PHYSICAL  (z)  PLANE 


FIG.2  MAP  IN  HOOOGRAPH  (w)  PLANE 
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FIG.  3  MAP  IN  W  PLANE 


FIG. 4  MAP  IN  £  PLANE 
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FIG.  5  INCOMING  JET  WITH  WALL 
ON  HIGH  PRESSURE  SIDE 


FIG.  6  PARTIALLY  CHANNELLED 
INCOMING  JET 
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